function y = Cauchy_CDF(x,m,b)
y = 0.5 + (1 / pi) .* atan((x - m)./b);
end